Pascal Principle for Diffusion-Controlled Trapping Reactions. 
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In this paper we analyse the long-time behavior of the survival probability Pa(^) of an A 
particle, which performs lattice random walk in the presence of randomly moving traps 
B. We show that for both perfect and imperfect trapping reactions, for arbitrary spatial 
dimension d and for a rather general class of random walks, PA(t) is less than or equal to 
the survival probability of an immobile target A in the presence of randomly moving traps. 
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Blaise Pascal has once asserted that all misfortune of 
man comes from the fact that he does not stay peace- 
fully in his room [1]. Taking this statement out of its 
philosophical context, it might be tempting to evoke it 
as the "Pascal principle" in regard to the problem of sur- 
vival of an A particle, which moves randomly in a sea 
of randomly moving traps B, (presented at mean den- 
sity "6"), and is annihilated upon the first encounter 
with any of them. This complex problem, which is inti- 
mately related to many fundamental problems of statisti- 
cal physics, has been attracting a great deal of attention 
within the last three decades (see, e.g., Refs. [2-4] and 
references therein). Many important results have been 
obtained but an exact solution is lacking as yet. 

In a recent paper [5], which focused on the behavior 
in the particular case when both A and Bs perform con- 
ventional diffusive motion, it has been claimed that it is 
intuitively clear that when the traps are initially (statisti- 
cally) symmetrically placed with respect to an A particle, 
the latter will on average survive longer if it stays still 
than if it diffuses. This represents, if true, a clear illus- 
tration of the " Pascal principle" , and implies that the A 
particle survival probability PA(t) obeys the inequality 

P A (t)<P A (t), (1) 

where P A (t) stands for the corresponding survival prob- 
ability of an immobile target A in presence of diffusive 
traps. The latter can be evaluated exactly [2,6]. Further- 
more, following the reasonings of earlier works [2,7], the 
authors of Ref. [5] constructed a lower bound on Pa(^) 
and showed that these bounds converge as t — > oo for 
systems of spatial dimension d < 2, defining thus the 
large-i asymptotic form of Pa(*) exactly! Subsequently, 
in Ref. [9], the arguments of Ref. [5] have been gener- 
alized to arbitrary symmetric random motion with par- 
ticles' and traps' trajectories characterized by a fractal 
dimension d u (not necessarily equal to 2, which is spe- 
cific for conventional diffusive motion), and here exact 



large-i asymptotic forms of Pa(*) have been obtained for 
systems of spatial dimension d < d w . 

We note that the inequality in Eq.(l) has been derived 
previously for the process of hopping transport of an ex- 
citation on a disordered array of immobile donor cen- 
ters in presence of randomly placed, immobile quenchers 
[10]. On the other hand, Eq.(l) is compatible with re- 
cent results on ballistic A + A — > annihilation process 
[11]. For trapping A + B — > B reactions (TR) involv- 
ing diffusive species the authors of Ref. [5] were unable, 
however, to prove the inequality in Eq.(l), but furnished 
only some arguments in favor of it. Consequently, results 
on exact asymptotic behavior of Pa (t) depend crucially 
on whether the inequality in Eq.(l) is indeed correct. 

In this paper we analyse, in the lattice formulation of 
the model, the Pascal principle- like inequality in Eq.(l). 
Following the line of argument of Refs. [6] and [8], we 
show that for both instantaneous (perfect) and imperfect 
TRs, for arbitrary spatial dimension d and for a rather 
general class of random walks (not necessarily conven- 
tional diffusion), P A {t) of a mobile A particle in presence 
of mobile traps, is less than the survival probability of 
an immobile target A in the presence of mobile Bs. Our 
proof assumes essentially that Bs perform a space- and 
time-homogeneous, unbiased jump process on the lattice 
sites. Behavior in the continuous-space systems, which 
requires a somewhat more delicate analysis (especially 
for d > 2), will be presented elsewhere [12]. 

Consider a d-dimensional hypercubic lattice contain- 
ing M sites. A single A particle is initially located at 
the origin, while N traps B are placed on the lattice 
at positions Yq ^ 0, where the superscript i here and 
henceforth numerates the traps, i = 1, . . . , N. 

In regard to particles dynamics, we suppose that the 
A particle performs a continuous-time jump process on 
the lattice sites and that the time interval between the 
consecutive jumps is a random variable. We denote 
then as F^ the A particle trajectory recorded at inte- 



1 



ger time moments k — 0, . . . , n, such that Ta = {Xo 

II. X i X„}, where is the vector of a lattice site at 

which the A particle resides at time moment k. Note that 
since the time interval between the consecutive jumps is 
a random variable, two successive positions are not nec- 
essarily different, and not necessarily nearest neighbors. 

Now, we suppose that the Bs perform identical and 
independent discrete time random walks: that is, at each 
tick of the clock each B can jump with a given proba- 
bility from a lattice site Y to another (not necessarily 
neighboring) site Y', or it may also remain at Y. We 
define then as the trajectory of the i-th B particle, 

= {Y£\ YP, . . . , Y ( n ] }, where denotes the po- 
sition of the i-th trap at time moment k, k = 0, 1, . . . , n. 

Next, let P(Y^\Y^ ] ) be the conditional probability 
of finding the i-th trap B at site Y^ at time moment 

(i) 

n, knowing that it started its random walk at Y . We 
assume now that random walks executed by the B par- 
ticles satisfy the following, quite general conditions: 

(i) the random walk is space- and time-homogeneous: 

P(Y«|Y«)=P(Y«-Y«|0), (2) 

(ii) at any time moment n and for any i [13]: 

P(Y« + Y«|Y«) < P(Y« = YW|Y«) = R n , (3) 
i.e., the conditional probability P(Yi J) ^ Y^|Y^) of 

(i) 

finding at time moment n the i-th trap at site Y„ dif- 
ferent from its starting point Yq ^ , is less than or equal to 
the probability P(Y ( r ? = Y^Y^) of finding it at time 
moment n exactly at the starting point; here, R n denotes 
the return probability of the random walk executed by 
the traps. By convention, Ro = 1. 

Finally, we consider two situations with respect to re- 
action; namely, when a) the A particle gets annihilated 
with probability p = 1 upon the first encounter with any 
of Bs (perfect reaction) and b) when the annihilation of 
the A particle takes place with probability p < 1 when 
an A and any of Bs occur at the same site (imperfect 
reaction). For computational convenience, we stipulate 
that for both situations reaction can take place only at 
integer time moments; that is, if at a non integer time A 
jumps on a site which is occupied by any B, A survives 
till the departure of this particle or an arrival of another 
B. The probability that both A and B particles jump on 
the same site simultaneously is clearly equal to zero. 

Let Pn^r^JYg 1 '') denote the conditional probability 
that for a given realization of the A particle trajectory 
T^, the i-th B particle starting its walk from the point 
Yq*' does not destroy (encounter) A up to time n. Since 
Bs move and act independently of each other, the condi- 
tional probability ^(r^KY^}) that, for a given real- 
ization of the A particle trajectory Ta and a given set of 



the starting points {Yq ^}, the A particle survives up to 
time n, is determined by 

^ n (T A \{Y^}) = f[P^(T A \Y^), (4) 

i=l 

and hence, the A particle survival probability obeys 

Pa{u) = (h n (T A \{Y^})) ) , (5) 

the average being taken first over the starting points of 
Bs and then over all possible trajectories T a- 

Now, supposing that Bs were initially uniformly dis- 
tributed on the lattice (excluding the origin) and droping 
the superscript "i", one finds from Eqs.(4) and (5) that 



P A {n) = 




Turning next to the thermodynamic limit, i.e. setting 
N, M = oo with a fixed ratio b = N/M, we find 




Consequently, the survival probability Pa(ti) can be 
thought of as the generating function of the probability 



Q„(r A |Y ) = l-P„(r A |Y ) (8) 

that for a given Ta, a single B, being at Y ^ at n = 0, 
destroys the A at some time moment < n. Note also that 
when the A particle is immobile, Eq.(7) reduces to 

P' A {n) = ex P J -b ^ (l-P„(0|Y )) 1, (9) 

[ Y #0 J 

which can be evaluated explicitly [2,6]. 

We seek now an upper bound on the survival proba- 
bility in Eq.(7). Let -Ffc(r,i|Yo) denote the conditional 
probability that a single B particle, being at Y at k = 0, 
encounters A for the first time at time moment k, given 
the A particle trajectory Ta is fixed. Then, the condi- 
tional probability Q n (TA\Y ) that a single B particle, 
starting from Y , destroys A at or before k = n obeys 

Q n (T A \Y )= F kFA\Y ). (10) 

0<fe<n 

Now, the conditional probability that the trajectory of B 
(extended after the possible annihilation of A) meets T A 
at time n (not necessarily for the first time) is clearly 
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P(Y n =X„|Y )-F„(r A |Y ) + 

+ Y P(y n = X n \Y k )F k (T A \Y ). (11) 

0<fc<n 

Summing both sides of the last equation over all initial 
positions Y ^ we obtain 

Y P(Y n =X n \Y ) = K n (T A ) + 

Y o #0 

+ Y P(Y n = X n \Y k )K k (T A ), (12) 

0<fc<n 

where K n (T A ) is defined as 

K n {T A ) = Y, F npA\Y ). (13) 

Yo/0 

We note that K n (T A ) has a meaning of a time-dependent 
reaction rate; using Eq.(13) we rewrite Eq.(7) as 



P A {n) = 




On the other hand, the survival probability P' A {n) of an 
immobile A particle, Eq.(9), can be written as 

P A (n) = cxpl-b Y K k \ =exp{ -6(5 n -l)}, 

[ 0<fe<n J 

(15) 

where K k = K k (T A = 0), while S n is the expected num- 
ber of distinct sites visited by a single B up to time mo- 
ment n (see, e.g., Ref [14] for more details). The last 
quantity is obtained directly by inversion of its generat- 
ing function S = ^2^Lq S n ^ n , which can be evaluated 

explicitly, S = (1 - £)~ 2 -R~S R = En=o R n£ n being the 
generating function of the return probability R n . 

We turn now to the comparison of K k (T A ) and K k . 
Using the normalization J2y P(Y n = X„|Y ) = 1, and 
the condition (i), we have Ey #o P(Y n = X„|Y ) = 
1 — P(Y n = X„|0). Consequently, in virtue of (ii): 

Y P(Y n =X n \Y )>l-R n . (16) 

Y #0 

Now, from Eqs.(12) and (16) we get 

K n (T A )+ Y P(Y n = X n \Y k )K k (T A )>l-R n . 

0<k<n 

(17) 

On the other hand, the inequality in Eq.(3) implies that 
P(Y„ = X„|Y fe )< J R n _ fe . (18) 



Recollecting that Rq = 1 and making use of Eq.(18), we 
thus enhance the inequality in Eq.(17), which now reads 

Y R n -kK k (T A )>l-R n , (19) 

0<fc<n 

and also becomes an equality when A is immobile. 

Further on, multiplying both sides of Eq.(19) by £™ 
and performing summation over n, we get 

K(Y A ) = Y K n (T A )C > — -g - I- (20) 

n=0 I 1 - Z,)R 

Next, taking into account that 

s(t a ) = yc( Y =t^t> ( 21 ) 

n=0 yo<fc<ri J ? 

we find from Eq.(20) the following inequality: 

(22) 

which implies that the generating function of the expres- 
sion in the exponent in Eq.(14), describing the A particle 
survival probability in case when it " leaves the room" and 
changes its position with time, is always greater than or 
equal to the generating function of the expression in the 
exponent in Eq.(15), which applies to the case when the 
A particle stays peacefully at its initial position. 

Hence, turning to the limit £ — > 1~ (n — > oo) and 
making use of the Tauberian theorems [15], we arrive at 
the conclusion that the desired inequality in Eq.(l) holds 
in the limit n — > oo. As a matter of fact, it can be 
shown that this inequality holds generally for arbitrary 
finite n; the proof in this statement is, however, rather 
cumbersome and will be presented elsewhere [12]. 

Lastly, we briefly outline the steps involved in the 
derivation of Eq.(l) in the general case when reaction be- 
tween an A and any of Bs is not instantaneous, but takes 
place with some finite probability p. Following Ref. [14], 
we suppose that here each trap bears "a gate", which 
may be either open or closed; in the former case the trap 
is reactive and annihilates the A particle upon the en- 
counter, while in the latter case it is inert with respect 
to reaction. The state of the gate on the i-th trap is 
characterized by a random variable Q such that d = 1 
(open gate) with probability p, and Q = (closed gate) 
with the probability 1 — p, respectively. Each Q updates 
its state at each tick of the clock; the updating process 
proceeds completely at random, without memory in time 
and without correlations with the gates imposed on other 
B particles. As shown in Ref. [14], such a model with 
stochastic, two-state gates corresponds to situations in 
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which the elementary reaction act is characterized by a 
finite intrinsic reaction constant K e i = p/(l —p). 

Now, we notice that Pa{ti) in this case can be still 
written in the form of Eq.(7) with Q n (TA\Yo) defined 
by Eqs.(8) and (10), but here F k (T A \Y ) = i> } (r A |Y ) 
should be interpreted as the conditional probability that 
the B particle encounters the A particle for the first time 
at time moment k exactly and moreover, that at this 
moment of time the B particle is in reactive state; the 
superscript "(p)" will signify that here we deal with im- 
perfect TR. Further on, let P (p) (Y n = X„|Y ) be the 
conditional probability that the trajectory of B meets 
r^4 at time n (not necessarily for the first time) and at 
this time moment B is in reactive state. For the model 
under study, such a probability obeys 



P ip \Y n = X„|Y ) = p P(Y n = X„|Y ) 



(23) 



and 



P W (Y™ = X n |Y ) = F„(r A |Y ) + 

+ ^ (p) (Y„-X„|Y fe )^(r A |Y ). (24) 



0<fc<n 



Summing both sides of Eq.(24) over Y ^ 0, we get 

P (p) (Y„-X„|Y o )-0)>p(l- J R„), (25) 

Y o #0 

and hence, we find 

P(l-Rn) <K<*)(T A ) + p Rn-k4 P) (T A ), (26) 



0<fe<n 



where Ki p) = ^ Yo#0 F^>(T A \Y ). Multiplying both 
sides of Eq.(26) by £™ and then summing it over n, n — 
0, . . . , oo, we arrive at the following inequality: 



S {P) (T ) > 



(I-?) 2 (l-p + pR) 



(27) 



where S^{T A ) is the generating function of the sum 
Eo< k < n K^(TA) = b-Hn(l/P A p \n)), P A p) (n) being 
the A particle survival probability for imperfect TR. 

Turning to the limit £ — > 1~ (n — ► oo), we notice that 
here (1 — £)i? <C 1 and hence, in this limit only the first 
term in the square brackets matters. On the other hand, 
this leading term S {p) = p(l - £)~ 2 / (l - p + pRj co- 
incides exactly with the expression obtained earlier [14] 
for the generating function of the exponent of the sur- 
vival probability of an immobile target A in presence of 
stochastically gated traps. The Tauberian theorem [15] 
then insures that also in this general case of imperfect 
TR the inequality in Eq.(l) holds as n — > oo. 



To conclude, we have proven here that in the long-time 
limit the survival probability of an A particle performing 
random walk on the sites of a d-dimensional lattice in 
presence of randomly moving traps is less than or equal 
to the survival probability of an immobile A particle in 
presence of randomly moving traps. This result holds 
for quite a general class of random walks as well as for 
perfect and imperfect trapping reactions. 

The authors thank J. Piasecki for fruitful discussions 
and for pointing us on the Pascal's assertion. We also 
acknowledge discussions with A. J. Bray and R.A.Blythe. 
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